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Abstract 

Alkali halide (100) surfaces are anomalously poorly wetted by their own melt at the 
triple point. We carried out simulations for NaCl(lOO) within a simple (BMHFT) 
model potential. Calculations of the solid-vapor, solid-liquid and liquid-vapor free 
energies showed that solid NaCl(lOO) is a nonmelting surface, and that the in- 
complete wetting can be traced to the conspiracy of three factors: surface anhar- 
monicities stabilizing the solid surface; a large density jump causing bad liquid-solid 
adhesion; incipient NaCl molecular correlations destabilizing the liquid surface, re- 
ducing in particular its entropy much below that of solid NaCl(lOO). Presently, we 
are making use of the nonmelting properties of this surface to conduct case study 
simulations of hard tips sliding on a hot stable crystal surface. Preliminary results 
reveal novel phenomena whose applicability is likely of greater generality. 
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1 Introduction 



Most liquids in nature are expected to wet very well the surface of their own 
solid at the melting point; and in fact they generally do. However, that is not 
always true and some examples of nonwetting do exist [1]. A striking case in 
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point is sodium chloride, where the NaCl(lOO) solid sTirface is badly wetted by 
its own melt, with an anomalously large partial contact angle of 48°(Fig. 1) [2]. 
How could one have predicted that, and what are the physical consequences? 
Moreover what applications, and what new phenomena, for instance in the 
field of friction, could this kind of non-self-wetting lead to? We addressed 
these questions using simple theory - basically surface thermodynamics and 
applied statistical mechanics - and computer simulation, in particular classical 
molecular dynamics (MD). 




Consider a NaCl droplet deposited on NaCl(lOO). From Fig. 1, 7sv = 7sl cos ^sl+ 
7lv cos ^lv and there will be partial wetting by the droplet if the three in- 
terface free energies satisfy 7sv < 7sl + 7lv (^sl — in our case). To ob- 
tain a microscopic theory, we must be able to calculate these three interface 
free energies. For the liquid surface, there is a well established method, the 
Kirkwood-Buff virial formula [3,4], that builds the surface stress (surface ten- 
sion), equal to the interface free energy, from mechanical variables, i.e., from 
forces and coordinates that can in turn be extracted from MD simulations. 
Nothing so straightforward is available for the solid- vapor and for the solid- 
liquid interfaces. The problem is that the solid can support stress, so that 
the surface stress, which is mechanical and easily calculable, and the surface 
free energy are no longer the same. For the solid-vapor interface, the prob- 
lem is solved by resorting to two parallel thermodynamic integrations carried 
out for two systems, identical in every respect, i.e., same number of parti- 
cles, same crystal structure, same cell etc., except that one system is a bulk 
without surfaces, the other is a slab with two (identical) surfaces. Integration 
of ///^^(T')rf(^) from some low temperature Tj up to T gives the free en- 
ergies of the two systems. The surface free energy is obtained at each T as 
half their difference (the slab has two surfaces) . For the sohd-hquid interface 
finally, not even that works, because the interface is stable just at T = Tm. 
Here, we circumvented the difficulty by turning so to speak experimentalists. 
We simulated a NaCl nanodroplet onto NaCl(lOO), let it settle, and extracted 
the contact angle ^lv [5]- Once ^lv, 7sv and 7lv are known, then 7sl can 
easily be obtained from Young's force balance equation given above. 
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Fig. 2. Calculated NaCl surface free energies, a) Solid-vapor 7sv- Note the metasta- 
bility up to about 150 K above T^. Dots: effective harmonic approximation, b) 
Liquid-vapor free energy: experimental and simulated values. 



Calculations and Results 



All results were produced using classical MD simulations, mostly in slab ge- 
ometries. Na"*" and Cl~ ions were assumed to interact via BMHFT poten- 
tials [6]. These potentials are standard and well tested, but they have a long 
range Coulomb part. This reduces dramatically the sizes of systems we can 
simulate - we used between 1500 and 5000 NaCl molecular units - and the 
simulation times we could afford. Luckily, two hundred psec were usually more 
than enough to attain equilibrium configurations at our high temperatures. 
Wc ran a bulk simulation first, and extracted the bulk melting - the tem- 
perature where the BMHFT solid and liquid coexist - at 1066 K, extremely 
close to the experimental 1074 K of NaCl [3,4]. 

We then simulated solid NaCl(lOO), and found that it remained dry and crys- 
talline up to 1210 K, well above the bulk T^. This kind of nonmelting behavior, 
where the solid surface can survive in a metastable state even above Tm, is 
exactly what can be expected when the hquid does not wet the solid [1]. The 
same slab system was then melted to generate a liquid slab - a liquid with 
two liquid- vapor interfaces. From these simulations, we extracted as explained 
earlier 7sv and 7lv shown in Fig.2. Here came three surprises. The first was 
the large decrease of 7sv with increasing T, strongly stabilizing the solid sur- 
face at high temperature. This is in part an anharmonic effect, and stems 
from the exceptional stability of long range Coulomb systems against large 
vibration amplitudes - even in bulk alkali halides largely exceeding the Lin- 
demann melting criterion [4]. The second surprise was to find 7lv so high - 
numerically equal to 7sv at Tm, when usually the solid surface is energetically 
much more expensive than the liquid one. The third surprise finally was that 
the surface entropy —dj/dT is a factor nearly three lower in the liquid than 
in the solid. This is strikingly contrary to our perception of a very disordered 
liquid surface (Fig. 3) as opposed to an ordered solid surface. As it turns out 
these two last surprises are related to one another, and are related to incipient 
molecular NaCl pair bond formation in the outer layer of the liquid surface. A 
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Fig. 3. Snapshot of simulated liquid -vapor interface at 




Fig. 4. Initial (a) and final (b) snapshots of NaCl droplet deposition on NaCl(lOO) 
at r = Tin. For the last 100 ps the drop remains in a metastable state forming a 
negligible internal contact angle ^SL) and an external contact ^lv angle of 50°. 

way to assess quantitatively this effect was found by recalculating the surface 
tension after removing molecular correlations, and showing that when this 
is done both 7lv (see ■jly in Fig. 2(b)) and the liquid surface entropy drop 
hugely [3,4]. 

Finally, we melted a small solid NaCl cube to form a nanodroplet, and de- 
posited it on the solid NaCl(lOO) slab surface (Fig. 4). From this final simula- 
tion we extracted an external wetting angle 6'lv of 50±5°, in good agreement 
with the experimental 48°. When inserted in the force balance equation along 
with 9sL ~ 0, that finally yields 7sl = 36 mJ/m^, which is about 1/3 of 7sv 
and 7lv. This relatively large value clearly demonstrates the poor adhesion of 
liquid NaCl to the solid, in turn reflecting the fact that they are very different - 
for example the liquid is 27% less dense than the solid. Nucleation studies [7,8] 
suggested even larger values of 7sl as large as 80 mJ/m^ at T ~800K. 

In conclusion, three separate pieces of physics conspire to cause the poor 
wetting properties of NaCl(lOO) by its own melt. The first is the extreme 
anharmonic stability of the solid surface itself. The second is the poor adhesion 
of the liquid to the solid. The third is the high liquid surface tension caused 
by the unexpected but very real surface entropy deficit [3,4]. 
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Fig. 5. Nanofriction simulation of a hard diamond tip ploughing into NaCl(lOO) near 
Tjn. The substrate is solid, but despite that the furrow heals away and recrystallizes 
in a very short time of a few ps. 

3 Nanofriction 



The next class of problems we are presently addressing, is as important as it 
is unexplored, is the high temperature nanofriction of a hard sliding tip near 
Tjn for a nonmelting surface like NaCl(lOO). On a normal solid surface - one 
that is wetted by its own liquid - a thin liquid film will nucleate below Tm 
causing a jump to contact with the approaching tip [9], and generally ruining 
it. A nonwetting surface like NaCl(lOO) will not do such thing. It is thus ideal, 
at least theoretically, to explore friction close to the melting point. Friction 
has been the subject of important simulation work [10,11]. High temperature 
friction is only recently beginning to receive some attention [12] and there is a 
clear need for more work. We intend precisely to explore and motivate future 
work in this extreme regime. 

We conducted preliminary sliding friction simulations for hard tips on NaCl(lOO), 
exploring both the deep ploughing regime with a sharp tip (Fig. 5), and the 
gentle grazing regime with a fiat tip. The tip was modeled to reproduce the 
interaction between rigid diamond and the NaCl surface[13,14]. For deep wear 
simulations we created a conical tip, composed of ~ 400 atoms. The diame- 
ter was around ISA and the height ~ 26A. On the contrary for light grazing 
wearless friction we use a completely blunt fiat tip composed of ~ 200 atoms 
(not shown). The fiat surface of the tip has a (111) plane of 13A diameter 
roughly. In ploughing regime the tip was moving with a constant speed of 50 
m/sec and the constant penetration depth of gA. The wear frictional force of 
Fig. 6 shows that at high temperature there is a regime where friction drops, 
although the surface is still fully solid. The tip scratches the solid surface, leav- 
ing behind a furrow which in this regime heals out closing spontaneously on a 
very short time of 10 ps. Direct inspection shows that the tip is surrounded in 
this high temperature regime by a seemingly liquid cloud, moving along with 
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Fig. 6. The dependence of friction force on temperature shows a decrease of nanofric- 
tion on NaCl(lOO) close to the melting point, as in ice skating. 



it and lubricating its motion. In this way, the tip effectively skates over the 
solid NaCl(lOO) surface. Discussions of skating near the melting point appear 
to exist so far only for ice [10]. 

Preliminary results show that under grazing conditions, the high temperature 
behavior of sliding friction is just the opposite [14]. The initially very low 
sliding friction of the flat tip over the cold surface rises with temperature, 
culminating in a frictional peak very close to the melting point. Here it appears 
that the large compliance of the (nearly unstable) hot NaCl surface lattice 
is responsible for the increased high temperature friction. There is to our 
knowledge no prior established mechanism of this kind for high temperature 
frictional increase. 



These predicted nanofriction phenomena are currently under closer scrutiny 
for further physical characterization. It is hoped that experimental efforts 
could be started to investigate their existence, which we believe to be more 
general for all nonmelting surfaces, rather than specific to the NaCl(lOO) sys- 
tem studied here. A more detailed report of this part of the work is forthcom- 
ing [14]. 
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